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Background: Nuclear effects can have a significant impact on neutrino-nucleus interactions. In particular, data 
from neutrino experiments with broad energy distributions require complex theoretical models that are able to 
take all the relevant channels into account as well as incorporate nuclear effects in both initial and final-state 
interactions. 

Purpose: We investigate neutrino and antineutrino scattering on iron and carbon in the energy range from 1 to 
50 GeV, which is relevant to current and coming experiments (MINOS, NOi/A, and Minerva). 

Method: The Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model, which implements all reaction channels 
relevant for neutrino energies under consideration, is used for an investigation of neutrino-nucleus reactions. 
Results: Our calculations are compared with the recent NOMAD and MINOS data for the integrated inclusive 
cross sections. Predictions are made for the differential cross sections for semi-inclusive final states (pions, kaons, 
and nucleons) for the MINOS and NOvA fluxes. 

Conclusions: Nuclear effects in the initial-state interactions may slightly change the inclusive nuclear cross 
section as compared to the free nucleon ones. Final-state interactions noticeably change the spectra of the 
outgoing hadrons. In the Minerva and NO^A experiments these effects should be visible in the kinetic energy 
distributions of the final pions, kaons and nucleons. Secondary interactions play an important role for strangeness 
production. They make it very difficult to extract the neutrino-induced strangeness production cross section in 
experiments using nuclear targets. 

PACS numbers: 25.30.Pt,13.15.+g,24.10.Lx 



I. INTRODUCTION 

Neutrino and antineutrino scattering on nuclei for neu- 
trino energies above 30 GeV was studied in several exper- 
iments starting in the 1980s (see Refs. for reviews). 
Based on muon detection, integrated cross sections were 
measured with overall precision of 2%; they grow lin- 
early with energy, which agrees with the predictions of 
the quark parton model. Double differential cross sec- 
tions with respect to muon variables were also measured, 
allowing one to extract the nucleon structure functions. 

Theoretically neutrino reactions at lower energies (a 
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few GeV to a few 10 GeV) still present a challenge be- 
cause here the methods of perturbative QCD are not re- 
liable and because contributions from quasi-elastic (QE) 
scattering, resonance (RES) production, background pro- 
cesses, and deep inelastic scattering (DIS) overlap. This 
is especially important in view of the broad energy dis- 
tributions of neutrino fluxes. The difficulty lies in choos- 
ing an appropriate way to describe the onset of the DIS 
processes, the so-called shallow inelastic scattering (SIS). 
This requires complex approaches that take all of the rel- 
evant channels into account. The lack of data at inter- 
mediate energies has long been an obstacle for a serious 
test of such approaches. 

New and coming experimental results are changing 
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the situation. Recently, the NOMAD collaboration per- 
formed measurements on an a composite target and re- 
ported the inclusive neutrino cross section on an isoscalar 
target for E v > 4.6 GeV with an accuracy of at least 4% 
Q . The MINOS experiment reported both neutrino (for 
E v > 3.48 GeV) and antineutrino (for E v > 6.07 GeV) 
cross sections also on an iron target with a compara- 
ble precision [J]. Cross sections at lower energies on 
(mostly) a carbon target will be measured by the NO^A 
experiment. The Minerva experiment intends to perform 
measurements on plastic (CH), iron, lead, carbon, water 
and liquid helium targets, which would directly allow to 
compare nuclear effects on various nuclei. In addition to 
muon detection, this experiment will also be able to re- 
solve various final states by identifying the tracks of the 
outgoing hadrons. 



On the theory side only very few papers have tried to 
describe not only the very high, but also the intermediate 
energy region. One of the most complete ones is the 
work by Kulagin and Petti [jj (see also Haider et al. Q]), 
who have presented a detailed description of neutrino 
inelastic inclusive scattering in terms of nuclear structure 
functions. In their work these authors also pay attention 
to a number of nuclear effects, such as Fermi motion and 
binding, and how these affect the extraction of structure 
functions. Our work is similar in spirit, but we aim for 
a practical implementation in an event generator that 
allows one to calculate not only inclusive cross sections, 
but also semi-inclusive ones as they will be measured by 
Minerva, for example. 



In this paper we study neutrino and antineutrino scat- 
tering on iron and carbon. Our results are compared with 
the recent MINOS and NOMAD data for inclusive cross 
sections. Predictions are also made for the spectra of the 
outgoing hadrons for the MINOS and NO^A neutrino 
fluxes. 



II. GIBUU TRANSPORT MODEL 

The GiBUU model has been developed as a transport 
model for nucleon-, nucleus-, pion-, and electron- induced 
collisions from some MeV up to tens of GeV. Several 
years ago neutrino-induced interactions were also imple- 
mented for the energies up to a about 2 GeV 0, Q] and, 
recently, the GiBUU code was extended to describe also 
the DIS processes for neutrino reactions. 

Thus, with GiBUU it is possible to study various el- 
ementary reactions on nuclei within a unified frame- 
work [9]. This is particularly important for broad- 
beam neutrino experiments, which inherently average 
over many different reaction mechanisms. Relevant for 
the present investigation is also the fact that the method 
and code have been widely tested for photon-induced 
as well as for electron-induced reactions in the en ergy 
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regime from a few hundred MeV to 200 GeV 
]. To be applicable for such a wide energy range the 
GiBUU model uses relativistic kinematics throughout as 
well as relativistic nuclear dynamics as expressed, e.g., in 
way how the mean-field potentials and cross sections are 
treated. Also the transport equation itself is covariant; 
for further details on these latter points, see Ref. Q. We 
stress that the calculations reported in the present paper 
are done without any fine tuning to the data discussed 
here with the default parameters as used in the general 
GiBUU framework. 

As already mentioned, broad band neutrino and an- 
tineutrino interactions with nucleons may result in sev- 
eral different channels. This is evident when looking at 
the energy distributions (see Fig. [1| of the two experi- 
ments, NOj/A and MINOS, that will be discussed later. 
The MINOS neutrino and antineutrino fluxes here are 
those used by the MINOS collaboration 



13] to deter- 



mine the total cross sections, that is, the neutrino and 
antineutrino parts of the low-energy NuMI beam in neu- 
trino mode. The same beam is used by the Minerva 
experiment. The NOi^A flux is the neutrino part of the 
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14 mrad off-axis medium-energy NuMI beam in neutrino 
mode 1141 . 
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FIG. 1. (Color online) Energy distributions for NO^A neu- 
trinos (solid red line) [uj], MINOS neutrinos (dashed green 
line), and MINOS antineutrinos [3] (short-dashed magenta 
line) normalized to unit area. The long tails of MINOS neu- 
trino (up to 50 GeV) and MINOS antineutrino (up to 38 GeV) 
fluxes are not shown. 

At hadronic invariant masses below the single- 
pion production threshold, W < 1.08 GeV, only the 
QE processes vn — > fj,~p and up — > [i + n are possible. 
At W > 1.08 GeV the single-pion production channel 
opens. The biggest contribution comes from production 
through the A [P33(1232)] resonance. With increasing 
W, excitations of resonances with higher masses, followed 
by their subsequent decays, give increasing contributions. 
For the prominent A resonance, the A — > Nir channel 
saturates the width at the 99% level. For higher res- 
onances, two and more pions, as well as other mesons, 
can be produced in the final state. In the GiBUU code, 
electroweak excitations of 13 resonances are accounted 
for, whose electromagnetic properties are taken from the 
Il6| . while their axial couplings are 
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assumed to be given by PCAC and dipole form factors. 
These processes as well as the non-resonant background, 
which gives a noticeable contribution in the same W re- 
gion, are discussed in Ref. jl7|. 



individual resonances become less and less distinguish- 
able. Thus, at high energies neutrino scattering needs 
to be described as a DIS process in terms of quark and 
gluon degrees of freedom. Within GiBUU this is done 
via the pythia event generator version 6.4.26 [181 ] ■ In it 
cross sections are obtained using the descriptions for the 
hard partonic processes, known from the Standard Model 
and QCD. Those are convoluted with the quark distri- 
bution functions, for which various known parametriza- 
tions can be used; the default is CTEQ 5L, a leading- 
order fit. The model separates elastic, single-diffractive, 
double-diffractive, and non-diffractive event topologies. 
The determination of final states is done with the help of 
the string fragmentation according to the LUND-String 
model or cluster collapses for low-mass configurations. 

Setting for the moment the difficulty aside to describe 
each of the possible neutrino channels (resonances, back- 
ground, DIS) separately, there appears a problem of pos- 
sible double counting. In the SIS region, i.e., the re- 
gion in between the resonances and DIS, the same phys- 
ical events can be considered as originating from decays 
of high-mass baryonic resonances or from DIS. In this 
region we switch smoothly from the resonance picture 
to a DIS description. While the resonance and single- 
pion background contributions are smoothly switched off 
above Wres-i the DIS contribution is switched on above 
Wdis-i; 



CTRES = CTres(M/) h/res _ 2 _ w/res _ 

for Wres-i < W < Wrbs-2, 



_ _ /TT/\ W-Wr>IS-l 

for Wdis-i < W < Wdis-2 



(1) 



High-mass resonances overlap: with increasing W, the 



The default parameters are W-res-i = 2.0 GeV, 
Wres-2 = 2.05 GeV, W mS -i = 1.6 GeV, and W ms _ 2 = 
1.65 GeV. The transition parameters given have been de- 
termined by comparison with the electroproduction data, 
as described in Ref. Q . The DIS processes in the interim 
region above Wdis— l and below Wres-2, where both 
resonance excitations and DIS are present, account for 
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the background processes giving a few mesons in the fi- 
nal state beyond the single-pion background. They also 
account for those (very few) resonances, which are not in- 
cluded in the GiBUU implementation because their elec- 
tromagnetic properties are not known. With this choice, 
the DIS events become noticeable at neutrino energies 
above about 1 — 2 GeV. At this point it should be men- 
tioned that in the absence of detailed data the cross sec- 
tions in the SIS energy regime carry some uncertainties; 
in particular, the 2tt channel, which clearly shows up in 
electron-induced inclusive cross sections, is completely 
undetermined, except for some old, low-statistics data 
13 . |20| . These uncertainties affect the total cross sec- 



tions on the 3% level. 

The GiBUU method to account for nuclear effects, 
which is described in detail in [9(, essentially factorizes 
the initial, primary interaction of the incoming neutrino 
with a target nucleon (or a pair of nucleons for 2p2h 
interactions) and the subsequent final-state interactions 
(FSI). The target nucleons are moving inside a mean-field 
nuclear potential with momenta determined by the Fermi 
motion in the local Thomas-Fermi approximation; all the 
potentials, defined in the local rest frame, have proper 
Lorentz-transformation properties and relativistic kine- 
matics are used throughout. While the relativistic Fermi 
gas (RFG) model has difficulties to describe details of 
the interactions of electrons with nuclei at lower energies 
and energy transfers in the QE region we have shown in 
Refs. [171 ? ] that at higher energies and energy transfers 
the GiBUU model leads to a description of QE scattering 
of the same quality as that reached in other more refined 
models of nuclear structure. We attribute this success 
partly to the use of a local RFG model as well as to 
the use of a momentum-dependent nucleon potential that 
takes nuclear binding effects and final-state interactions 
in the first, primary interaction into account. In the en- 
ergy region treated in this paper QE scattering amounts 
anyway only to a small part of the total cross section 
while pion production and DIS become more important. 



For these latter two processes GiBUU gives results that 
are very reliable when compared to data (see Figs. 28, 
29, and 33 in Ref. fl). 

Thus, we consider the cross sections of neu- 
trino/ antineutrino interactions 

v{k^)N(p>") -> fi-(k'^)X , 

with a bound nucleon, which has a non- vanishing three- 
momentum and effective mass m* N that takes into ac- 
count the nuclear potential. In the interaction possible 
Pauli blocking of the final state is explicitly taken into ac- 
count. The cross sections are then obtained by summing 
over the interactions with all nucleons taking the flux fac- 
tor (q ■ p) I '(qoPo) , which transforms the neutrino- nucleon 
flux into the neutrino-nucleus flux 
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into account. 



The total absorption cross section is then calculated 
as CT t ot = ctqe + ores + 0"bgr + cdis- For QE scattering 
and resonance excitations the primary interaction of a 
neutrino with an individual target nucleon is described 
taking into account the full in-medium kinematics with 
an energy-dependent nuclear potential for both the in- 
coming and the outgoing baryons. The background cross 
section is essentially the difference between the experi- 
mentally known single-pion production cross section and 
the resonance contribution to these processes. The de- 
tails are given in Ref. citeLeitner:2008ue. 

High-energy electron and neutrino data on nuclear tar- 
gets, which are dominated by the DIS processes, have 
conventionally been analyzed by assuming that the in- 
clusive double-differential cross sections assume the same 
form as those for the free reaction and all possible nu- 
clear effects have been absorbed into the nuclear modifi- 
cation of the structure functions (see, e.g., Refs. 
In contrast, the present model works quite differently. 
The calculations explicitly contain nuclear effects such as 
Fermi motion and Pauli blocking, but no intrinsic struc- 
tural change of bound nucleons. Any disagreement of 
our curves with the data will thus reveal "genuine" nu- 
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clear in-medium effects not included in the model. For 
QE scattering and resonance excitations the effects of the 
nuclear potential are explicitly taken into account by us- 
ing spectral functions with the properly shifted real parts. 
GiBUU also allows us to account for the in-medium width 
modification of the resonances; in the present high-energy 
calculations, however, we neglect this effect for simplicity 
and because its influence on the observables calculated 
here is small. 

pythia, used to generate the DIS events, is a free- 
nucleon generator, which does not allow explicit use of 
a nuclear potential. We try to account for this and thus 
determine the elementary cross section appropriate for a 
bound nucleon by adjusting the input kinematical vari- 
ables for PYTHIA. Various prescriptions to do this have 
been discussed in Ref. We will come back to this 
point later in Sec.|V] Here we just mention that the anal- 
ogous uncertainty also exists in the standard treatment 
of electron- and neutrino-induced inclusive cross sections 
where one has to decide at which kinematical variables u> 
(energy transfer) and Q 2 (four-momentum transfer) the 



structure functions have to be read off 
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After being produced in the initial interaction, out- 
going hadrons propagate throughout the nucleus. In 
GiBUU this process of FSI is modeled by solving 
the semi-classical Boltzmann-Uehling-Uhlenbeck equa- 
tion, which, again, respects relativity jgj]. It describes 
the dynamical evolution of the one-particle phase-space 
density for each particle species under the influence of 
the mean field potential, introduced in the description 
of the initial nucleus state. Equations for various parti- 
cle species are coupled through this mean field and also 
through the collision term. This term explicitly accounts 
for changes in the phase-space density caused by elastic 
and inelastic collisions between particles. We note here 
that, contrary to some other uses of this term, we call 
FSI all the interactions that take place after the initial, 
primary reaction while the hadrons propagate through 
the nucleus. However, the same potential that governs 



this later evolution is also present during this first reac- 
tion and affects its outcome, mainly due to the energy 
dependence of the mean-field potential that changes the 
final-state phase space. 

At higher energies and, in particular, higher Q 2 , for- 
mation times and color transparency phenomena may be- 
come important. Relevant data were taken by EMC and 
HERMES collaborations some years ago. In Ref. 12j we 



have analyzed these data in the energy regime from about 
10 to 200 GeV and have found that only a linear increase 
of the prehadronic cross sections with time fits both sets 
of data simultaneously. Therefore, the quantum-diffusion 
model of Farrar et al. 25| is implemented in GiBUU. At 



the lower energies formation times are determined by the 
lifetimes of resonances, which are explicitly propagated, 
and as such are automatically included in GiBUU. 

III. RESULTS ON THE NUCLEON 
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FIG. 2. (Color online) Cross section da/dQ 2 per nucleon for 
CC neutrino scattering off an isoscalar target for E v — 4 GeV. 

At low neutrino energies up to about 2 GeV, the cross 
sections for neutrino reactions on the nucleon were pre- 
sented in Ref. p| and showed good agreement with the 
experimental data. At higher energies (above about 
E v « 2 — 3 GeV) QE, resonance (A and higher reso- 
nances) and background (bgr) contributions practically 
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do not change with energy. This is because they strongly 
fall off with increasing Q 2 , as shown in Fig. [2] for the 
example of an isoscalar target for a neutrino energy of 
4 GeV. The DIS contribution, on the other hand, is 
much flatter in Q 2 ; its absolute value will also grow with 
increasing energy. 

Thus, at neutrino energies above a few GeV the major 
contribution to the total cross section comes from DIS. 
According to the predictions of the parton model, the 
DIS cross section grows linearly with energy. At high 
neutrino energies the data are, therefore, conveniently 
presented as cross section per energy a^ ^/E v . 

This cross section was measured starting in the 1980s 
by the CCFR, BEBC, CHARM, IHEP, NuTeV, CHO- 
RUS and NOMAD collaborations; see Ref. [26( for the 



plot with all data and references. The world average val- 
ues derived for E v > 30 GeV are given for an isoscalar 
target: a^ Qt /E v = 0.667 ± 0.014 • I0~ 38 cm 2 /GcV for 
neutrinos and ^ Qt /E p = 0.334 ± 0.008 • lO" 38 cm 2 /GeV 
for antineutrinos. 

Figure [3] shows the results of our calculations of the 
total cross section as well as the DIS and other contri- 
butions indicated in the figure for both neutrino and an- 
tineutrino scattering on an isoscalar target. The world 
average values and their error bands are indicated for 
comparison for E v > 30 GeV. The NOMAD experiment 
has recently performed measurements on composite tar- 
get with a measured composition of 52.43% protons and 
47.57% neutrons [jj. The measurements were then cor- 
rected for the nonisoscalarity, and the results are pre- 
sented as isoscalar cross section. The data points are 
shown in Fig. [3] as solid triangles. 

For neutrinos, the DIS contribution becomes larger 
than the A contribution at about 3 GeV, and at 5 GeV 
it is already about 60% of the total and reaches 95% at 
higher energies. The rest is to be attributed to other 
channels. For antineutrinos, the DIS contribution be- 
comes larger than the A channel at about 4 GeV; it is 
40% of the total at E s ~ 5 GeV and reaches 95% at 



higher neutrino energies. This clearly has implications 
for theoretical analyses of the cross section measurements 
in the MINOS and NOz/A experiments. 

The dip in the total antincutrino cross section at E v ~ 
3 — 4 GeV (there is an indication for a similar effect in 
the neutrino cross section) shows up as such only in the 
l/E v scaled cross section. It is caused by an interplay 
of the downfall of the resonance contributions and the 
rise of DIS. It could indicate that our model misses some 
strength here. This intermediate energy region would be 
most sensitive to the pion (one or more) background. A 
clarification must wait until much more precise data for 
l7r and 27T production become available. 

The figure also shows for both neutrinos and antineu- 
trinos a slight decrease of the calculated atot/E at large 
energies. In our calculations this slow decrease is due 
to the quasielastic and resonance contributions that fall 
with energy. The high-energy data also exhibit a small 
downward slope usually attributed to perturbative QCD 
and heavy quark effects. This latter slope was not taken 
into account in deriving the world average values, where 
it was assumed to be negligible [lj]. 

IV. QUASIELASTIC SCATTERING 

As we have shown in Fig. [3J with increasing neutrino 
energy the CC QE cross section gives smaller and smaller 
relative contribution to the total cross section. However, 
this channel is still of importance because of its special 
role in the reconstruction of the neutrino energy. The MI- 
NOS [27J, Miner j/a 28], and ArgoNeuT [29( experiments 



all aim at studying this channel. They use the same 
neutrino flux, but have different detectors and different 
criteria for selection of QE events, which, again, differ 
from those of low-energy experiments like MiniBooNE 
and K2K. 

Here we present in Fig. [4] only the QE cross section 
as a function of Q 2 , calculated for an axial mass of 
M A = lGeV. The plot is shown for the MINOS flux; 
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FIG. 3. (Color online) Total (red solid), DIS (blue dashed) and other cross sections per energy per nucleon for antineutrino 
(left) and neutrino (right) inclusive scattering on an isoscalar target. Curves are compared with the world averaged values, 
denoted by the dashed band. Data points from the NOMAD experiment [jj are shown as solid triangles. 
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FIG. 4. (Color online) Cross section da/AQ 2 per neutron 
for neutrino CC QE scattering as a function of squared four- 
momentum transfer. The cross section is averaged over the 
MINOS flux. 

the results obtained with the NO^A flux are nearly in- 
distinguishable, as expected for neutrino energies above 
1.5 GeV. The elementary QE cross section used here is 
given in Ref. [7| ■ The dashed line shows the free cross sec- 
tion on a single neutron; the solid line takes into account 
nuclear effects in iron. Those reveal themselves only at 
low momentum transfer; for Q 2 > 0.2 GeV 2 the "neu- 
tron" and " 56 Fe" curves are practically indistinguishable. 

We show the flux-integrated QE result here in order to 
illustrate that the Q 2 distribution for QE scattering looks 



very similar to that at the lower-energy experiment Mini- 
BooNE. Comparison with Fig. [5] shows that the different 
reaction mechanisms have quite different Q 2 behavior. 
Whereas both QE and the A contribution fall off quite 
quickly with Q 2 , driven by the Q 2 dependence of their 
form factors, the DIS contribution is, after an initial rise, 
rather flat out to larger momentum transfers. 

We stress here that the solid curve in Fig. 0] depicts 
the true CC QE cross section without any initial pion 
events and without any initial 2p-2h interactions. The 
latter would contribute not to QE scattering, but to the 
inclusive pionless cross section. Taking the MiniBooNE 
results for guidance they would contribute about 30% to 
the latter, but their overall influence on the total cross 
sections would be negligible at the higher energies where 
the dominating DIS cross section goes cx E v . 

V. NUCLEAR EFFECTS IN DIS 

All the neutrino experiments, used to derive the world 
average cross sections on the nucleon mentioned in the 
previous section, were actually performed on nuclear tar- 
gets. If the target was nonisoscalar, a corresponding cor- 
rection was introduced and the isoscalar cross section was 
extracted jlj]. Such a procedure as well as the introduc- 
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FIG. 5. (Color online) Cross sections rjtot/E v per nucleon for antineutrino £uFe — > fi + X (left) and neutrino v^Fe — > fi~X 
(right) inclusive scattering off iron. Experimental data are taken from Refs. [j] (MINOS, open circles) and (IHEP-JINR, 
solid diamonds). 



tion of the "world-average" value is meaningful only if 
nuclear corrections are very small. 

The EMC effect shows that there are nuclear correc- 
tions to the free cross sections for electrons [32[. For 



neutrino reactions, however, the situation is controver- 
sial. On one hand, nuclear parton distributions, based on 
electromagnetic scattering data and intended for descrip- 
tion of both charged lepton and neutrino reactions, were 
introduced. For a review and a list of recent parametriza- 
tion see, for example, Ref. 



recent investigation 



34, 



On the other hand, a 



35| showed that in neutrino re- 



actions nuclear corrections to parton distributions have 
about the same magnitude as for electrons, but have a 
very different dependence on the Bjorken-a; variable. The 
topic remains controversial, with the hope that future 
precise Minerva results on various targets will clarify the 
situation. 

As we have already mentioned above, the GiBUU code 
uses PYTHIA for the simulation of the DIS processes. 
Since the PYTHIA code was designed for elementary re- 
actions we have to provide some "quasi-free" kinematics 
as input to pythia that removes the effects of the bind- 
ing potential on the nucleon. Various prescriptions to 
do this have been used (for details see Ref. Q) and are 
compared with each other and with the free cross sec- 



tion in Fig. The corresponding cross sections are de- 
noted as "F-NO" (the invariant energy W 2 of the boson- 
nucleon system is calculated as (k — k' + p) 2 and not 
corrected), "F-CM" (bound nucleon is boosted into the 
center-of-momentum frame and the nuclear potential is 
removed from its energy, the nucleon then is boosted 
back, and W is calculated), and "F-MED" (W 2 is taken 



as (k - k' + p) 2 - 



1 *n 2 



m 2 N ). In all these calculations 



parton distributions appropriate for free nucleons have 
been used. 

It is seen that the difference between these various pre- 
scriptions is quite small (about 2 % at the lowest energy 
and less at the highest energy); also, the results all ap- 
proach the free cross section at the highest energy. We 
consider small differences between the results obtained 
with the various prescriptions mentioned above as intrin- 
sic uncertainty of the GiBUU code, reflecting the lack of a 
detailed understanding of nuclear effects. No other event 
generator, as far as we know, accounts for nuclear correc- 
tions in high-energy neutrino reactions. Nuclear parton 
distribution functions from Ref. 3(| are also implemented 
as one of the options to use. To avoid double counting, 
nuclear potential and Fermi motion are switched off in 
such calculations. The result ( "nuclearPDF" ) as well as 
the free cross section for iron composition ("26p+30n") 
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are also shown in Fig. [5] 
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FIG. 6. (Color online) Ratio of the antineutrino to neutino 
cross sections in scattering off iron. Also shown is the cal- 
culated ratio for an isoscalar-corrected target. Experimental 
data are taken from Ref. (MINOS, open circles). 

Figure [5] shows, that for antineutrinos our curves are 
within the spread of the MINOS and IHEP-JINR data 
31|. The overall agreement of our calculations with the 
data is, therefore, better than the agreement of the data 
with each other. At low energies, where the main con- 
tribution comes from the QE and resonance production, 
nuclear effects are known to reduce the neutrino and an- 
tineutrino cross section. This is why, at Ep < 5 GeV, 
the curves "F-NO", "F-CM", and "F-MED" that take 
into account the nuclear effects explicitly lie noticeably 
lower than the "26p+30n" curve. At high energies, how- 
ever, the curves converge toward each other. The "nu- 
clearPDF" curve, which implies modification of the DIS 
channel only, coincides with the "26p+30n" curve at low 
energies and consistently deviates from it to lower values 
of the cross section at higher energies where DIS domi- 
nates. The peak and dip in the region 3 — 4 GeV in the 
free cross section have the same origin as for the isoscalar 
cross section, as discussed in the previous section. Nu- 
clear effects, mainly Fermi motion, wash out this struc- 
ture so it is no longer visible in the nuclear cross sections. 

For neutrinos our curves are in good agreement with 
the recent MINOS experiment, they lie within the errors 



of the data points. For both neutrinos and antineutrinos 
nuclear effects are noticeable at low energies. For higher 
energies, E v > 5 GeV, the curves "F-NO" "F-CM" and 
"F-MED" all approach each other and the "26p+30n" 
curve, reflecting the expected disappearance of nuclear 
effects with increasing energy. 

In Fig. [5] we show the ratio of antineutrino cross section 
to the neutrino one as a function of energy, which is in 
good agreement with the recent MINOS data. The ratio 
rises with energy, gradually flattening out. It is below 
its asymptotic value obtained for an isoscalar-corrected 
target (« 0.5 at high energies |lj); the corresponding 
isoscalar curve is also shown in the figure. It is interesting 
that this ratio is remarkably insensitive to any nuclear 
effects, even at the lower energies, as illustrated by the 
fact that now all the curves for the various in-medium 
correction methods lie essentially on top of each other. 



VI. SEMI-INCLUSIVE CHANNELS 

In this section we present spectra for outgoing pions 
and nucleons. All calculations have been done for the 
MINOS flux on an iron target and the NO^A flux on a 
carbon target. As emphasized, for example, in Ref. 



37], 



acceptance cuts or detector thresholds can have a signifi- 
cant and non-trivial influence on the measured values. In 
present calculations, no such cuts and no detector thresh- 
olds are assumed for the outgoing hadrons. 

In many cases, FSI significantly modify the shapes 
of the final particle spectra. Particles get slowed down 
in the medium by collisions with other nucleons lead- 
ing to a pileup of cross section at low kinetic energies. 
Such modification is seen, for example, in photopion pro- 
duction ll| and is well described by GiBUU. A similar 
change should be observed in neutrino reactions. In the 
following we now discuss nucleon knock-out, pion pro- 
duction, and strangeness production in neutrino-induced 
reactions. 
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A. Nucleon kinetic energy distributions 

Knock-out nucleons may give an insight into the ac- 
tual primary interaction process and, in particular, into 
the question if the impulse approximation is sufficient or 
if multi-nucleon interactions are essential. The kinetic 
energy distributions of outgoing nucleons are shown in 
Fig. [7] for single-nucleon events (one nucleon of a given 
charge and no other nucleons in the final state) and in 
Fig. |S] for multi-nucleon events (at least one nucleon of 
the indicated charge and any number of other nucleons). 

The decrease at higher energies (T > 0.05 GeV) due to 
FSI is natural to expect, because a nucleon can rescatter 
in the nucleus and knock out another nucleon; the nu- 
cleon is then gone from the single-nucleon channel. At 
the same time, its kinetic energy would be spread be- 
tween the two secondary nucleons. If these two have 
an energy large enough, they could, in turn, produce 
more lower-energy nucleons. The same knock-out can 
be caused by a pion produced in a primary interaction. 
This process, which can develop as a cascade, leads to an 
increase of multi-nucleon events with low kinetic energies. 

This is indeed seen in Fig. [8] which shows the kinetic 
energy spectra for multi-nucleon knock-out for both pro- 
tons and neutrons. Here the cross sections before FSI 
are quite similar to those in the single-nucleon events, 
but there is a dramatic effect of FSI with a huge in- 
crease of the cross section at smaller energies (Tjv i$ 0.05 
GeV) . This steep rise at smaller energies was already ob- 
served earlier in Refs. 0, and is due to the energy 
loss connected with multiple scattering as just discussed. 
We note that this behavior can only be obtained in a 
transport-theoretical treatment of the neutrino-nucleus 
interaction; models that work with an optical potential 
for the outgoing nucleons can usually only describe their 
loss of flux, but not where this flux goes. An exception 
is a treatment based on Feshbach's projection formalism, 
that, in principle, is equivalent to a full coupled-channels 
treatment and, thus, allows also a description of inelastic 



processes and energy loss. Applications of this formalism 
have so far been restricted to inclusive cross sections [? ] 
and to exclusive proton knock-out reactions [? ]. 

For NO^A flux calculations, which are done for car- 
bon nucleus, the picture is similar. Since the nucleus is 
smaller, the hadrons can leave it earlier and the cascade 
processes have less time to develop. So, FSI lead to lower 
nucleon multiplicities (as compared to iron) and corre- 
spondingly higher kinetic energies. As a consequence, 
for single-nucleon events, the suppression is smaller than 
that for iron and for multi-nucleon events the rise at low 
energies is not so steep. 

At this point it is in order to comment briefly on the 
possible effects of initial many-body (so-called 2p-2h) ex- 
citations that are not contained in the calculations dis- 
cussed here. The presence of such excitations, which is 
well established in inclusive electron scattering in the re- 
gion between the QE peak and the first resonance (A) 
excitation, has recently been invoked to explain the dif- 
ference between the QE-like cross sections measured in 
the MiniBooNE experiment and theoretical calculations 



39-42]. Here we note that for the MINOS flux most of 
the multi-proton knock-outs come from DIS events and 
that already QE plays only a minor role. Any 2p-2h con- 
tributions thus necessarily would be quite small and can 
be neglected. 

B. Pion kinetic energy distributions 

Figure |H] shows the ir + , ir° and tt~ spectra for single- 
pion events for neutrino- and antineutrino-induced reac- 
tions. Single-pion events are defined as those having only 
one pion of a given charge and no other pions in the final 
state. 

For the dominant channels (n + production in neutrino 
reactions and ir~ in antineutrino ones), the FSI decrease 
the cross section at T n > 0.2 GeV. At lower neutrino 
energies this is mainly due to pion absorption through 
ttN A followed by AN -> NN. At higher ener- 
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and no other nucleons) production in neutrino and antincutrino scattering off iron and carbon. The left block of figures shows 
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FIG. 8. (Color online) Kinetic energy distributions per target nucleon for multi-nucleon (at least one nucleon of the indicated 
charge and any number of other nucleons) production in neutrino and antineutrino scattering off iron and carbon. 



gies pions can also be absorbed through rj and A pro- 
duction ttN — > rjA , production of higher resonances 
ttN -» R followed by RN NN or R -)• rjN, non- 
resonant pion absorption ttNN — > NN, production of uj 
mesons 7rA^ — > uiN, c/> mesons irN —> <pN, and strange 
mesons ttN — > EAT, AK, KKN. All these channels (and 



more) are included in GiBUU and contribute to pion ab- 
sorption. 

In addition, pion scattering in the FSI also decreases 
the pion energy. Here elastic scattering as well as DIS 
events of the type irN — > multi-TrA^ deplete the spectra 
at higher energies and accumulate strength at lower en- 
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FIG. 9. (Color online) Pion kinetic energy distributions per target nucleon for neutrino and antineutrino scattering off iron 
and carbon for single-pion production (one pion of the indicated charge and no other pions are produced). 



ergies. Thus, an increase of the cross sections is observed 
at T,,. < 0.15 GeV, where the cross sections after FSI are 
higher than before, and decrease above this energy. Ad- 
ditionally, low-energy pions may come from the reactions 
such as T) N — > R followed by R — > irN. Altogether this 
leads to a significant change of the shape of the spec- 
tra. In particular, there is a strong build-up of strength 
around = 0.06 GeV, where the cross section after FSI 
is about 50% higher than before. This is primarily due 
to the slowing down of pions by FSI and the low ir — N 
cross section in this energy region. We note that the size 
of this effect depends somewhat on the treatment of the 
collisional width of the A resonance j43j . 

Pions can also be produced in interactions of secondary 
nucleons. In addition, pion scattering can also involve 
pion charge exchange. For neutrino-induced reactions, 
the 7r + n — > tt°p scattering in the FSI is the main source of 
side-feeding into the 7r° channel, leading to a noticeable 



increase of the ir° cross section at low T ff . The same 
effect was found for low-energy neutrino reactions Q]. 
Now the cross section after FSI is about 200% higher 
than before. The inverse feeding is suppressed, because 
less 7T° than ir + are produced at the initial vertex. The 
same mechanism of side feeding from dominant to sub- 
dominant channel through ir~p —± ir°n is working also 
for antineutrino-induced reactions. 

For the least dominant channel (7r~ production in neu- 
trino reactions and n + in antineutrino ones) , the FSI (in 
particular, side feeding into this channel) represent the 
main source of the events observed. Indeed, the ir~ yield 
with neutrinos and n + yield with antineutrinos expected 
for the MINOS experiment (cf. Fig. |9]) is, at the peak, 
only a little lower than that for 7r° production. 

The corresponding results for the kinetic energy distri- 
bution of pions in multi-pion events, defined as consisting 
of at least one pion of a given charge and any number of 
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FIG. 10. (Color online) Pion kinetic energy distributions per target nucleon for neutrino and antineutrino scattering off iron 
and carbon for multi-pion production (at least one pion of the indicated charge and any number of pions of other charges are 
produced) . 



other pions, are shown in Fig.[TUJ The cross sections here 
are 2-3 times higher than those of single-pion production 
in Fig. |9l but the spectral distributions are very similar. 

The observed pions can be produced by different mech- 
anisms. It is, therefore, interesting to analyze where the 
various final states come from. Fig. Qj] shows the ori- 
gin of the pions (that is, the initial vertex at which the 
pion was produced) in the dominant channels for various 
final states. It is interesting to see that even for the MI- 
NOS flux, which peaks at 3 GeV and has a high-energy 
tail, single-pion production receives its major contribu- 
tion from A resonance production and its following de- 
cay. The second largest contribution comes from DIS. 
This reflects the fact that even in a DIS event the final 
state may involve a A that subsequently decays into a 
pion. For the other final states with more than one pion 
DIS dominates, but the A is still visible. This reflects the 



fact that the pion and nucleon produced in the decay of 
the primary A undergo FSI which result in several pions 
in the final state. The contribution from the QE vertex 
is very small but nonzero. In this case the outgoing pion 
can be produced only during the FSI, for example, due 
to the NN ->• NA scattering followed by A ->■ Nn. 

The integrated pion cross sections, averaged over the 
MINOS and NOz^A flux distributions, are shown at the 
end of this section in Table |U Notice here that numeri- 
cal integration of the single-pion distributions in Fig. [5] 
directly leads to the values presented in the table. For 
the multi-pion events this is, however, not so. Indeed, 
in the multi-pion distribution each pion in a given event 
gives a contribution to the results in Fig. [TO] Thus, for 
example, an event with 37r + in the final state is counted 3 
times in the multi-7r + distribution in Fig. 1101 For the in- 
tegrated multi-7r + cross section such an event is counted 
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FIG. 11. (Color online) Pion kinetic energy distributions per target nucleon for tt + production in neutrino and ir~ production 
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only once. 

C. Kaon kinetic energy distributions 

An advantage of the fine-grained detector used in the 
Minerva experiment is its ability to measure the energy 
loss in the calorimeter with high precision and, thus, to 
distinguish among pions, kaons and nucleons. Since the 
particle identification is based on time-of-flight measure- 



ments, the lower the kinetic energy, the better the sep- 
aration will be; the upper boundary for distinguishing 
kaons from pions is T « 500 MeV. As we have seen, for 
the low-energy NuMI beam we expect that the majority 
of the outgoing pions and nucleons lie below this energy 
region. A similar situation is expected for kaons. 

At neutrino energies near the kaon production thresh- 
old and slightly above, the cross section can be described 
in term of hadronic degrees of freedom and is expected to 
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be at the level of 10 -41 cm 2 for both kaon production in 
neutrino reactions and antikaon production in antincu- 
trino reactions 



44 



45]. In the experiments discussed 



here these strangeness-changing events are clearly out- 
numbered by kaons produced through DIS which domi- 
nates the kaon cross sections. Therefore, we have imple- 
mented only the strangeness conserving channels which 
are described by pythia in terms of quark and gluon de- 
grees of freedom as ss pair production with the following 
fragmentation. 

Before FSI, kaons can originate from events such as 
v + N — ► KA, vN — > K + K~N, or more inclusive events 
with some pions in the final state. Due to FSI a kaon 
can be produced in secondary pion rescattering on nu- 
cleons in processes such as irN — > Y,K, AK, KKN, and 
by nucleon rescattering via NN — > NUK, NAK. Once 
produced, kaons can also rescatter elastically, K + N — > 
K + N, and inelastically, K + N — > K + Nir, reducing their 
energy. In addition, they can undergo charge transfer, 
K + n — > K°p. Antikaons can only be produced together 
with kaons and they can be easily absorbed in reactions 
such as K~N Air~ or K~N — > £°7r~ which dom- 
inate the total K~N cross section at small momenta 
<« 0.8 GeV. Because these two processes proceed 
through anti-strange resonances these two cross sections 
become very large at small kaon momenta leading to a 
large absorption of slow antikaons. 

Kaons are interesting to study because K + have a long 
mean free path because of strangeness conservation; they 
can disappear only by charge exchange into K°. On the 
other hand, they could also be produced in secondary 
collisions. This makes them sensitive to possible forma- 
tion times of high-energy hadrons in the medium. In 
Ref. fiol ] it was shown that for kaon-photoproduction on 
nuclei, which is closely related to the neutrino-induced 
production investigated here, the final-state interactions 
actually increased the cross section for kaon production. 
In that same paper a significant dependence of the fi- 
nal kaons on the hadron formation times was also found. 



This was mainly due to secondary reactions of primary 
high-energy pions. Such pions have a very large Lorentz- 
boosted formation time in the nuclear target frame which 
suppresses their secondary interactions. As a conse- 
quence, for photon energies larger than about 3 GeV the 
kaon yield actually increased when the formation time 
for hadrons was set to zero. On the contrary, the an- 
tikaon yield decreased for zero formation time due to ear- 
lier absorption. A very similar result had been found in 
Refs. [12|:|46|]. We recall from the latter calculations that 
the so-called leading hadrons which carry quarks from 
the initial reaction partners are less strongly affected by 
the formation times than those hadrons that are being 
produced from the sea [12]. Among other effects this 
leads to a difference in the influence of FSI on kaons and 
antikaons. 

Fig. HU shows the kaon kinetic energy distribution for 
multi-kaon events (at least one kaon of a given charge and 
any number of other mesons). Comparison of Fig. [TU] 
and Fig. [T^] shows, that the ratio of kaons to pions is 
expected to be at the level of 0.03 — 0.1 (depending on 
kinetic energy) . This is very similar to the recent result of 
Na61/SHINE experiment on proton scattering off carbon 



47]. 



Since nuclei contain only positively charged and 
neutral nucleons, and any reaction should conserve 
strangeness, baryon number and charge, one would ex- 
pect noticeably fewer antikaons than kaons in neutrino- 
induced reactions even before FSI. The kinetic energy dis- 
tribution for antikaons (K~ and K°) is shown in Fig. [13] 
The corresponding cross sections are indeed 3-10 times 
smaller than those for kaon production in Fig. Q21 

Interesting here is that the small antikaon cross sec- 
tions do not show any significant reabsorption, contrary 
to our naive expectation. For the NO^A flux the antikaon 
output after FSI is nearly the same as that before FSI. A 
closer investigation of this result has shown that in this 
channel, where the cross section for the primary vertex is 
small, secondary production of antikaons is very impor- 



16 



v Fe -> K + 
multi-kaon production 




0.1 
0.08 
0.06 
0.04 
0.02 


0.08 
0.06 
0.04 
0.02 





vC^K + withF SJ, 
w/o FSI 



multi-kaon production 




0.5 1 



1.5 



2.5 



FIG. 12. (Color online) Kinetic energy distributions per target nucleon for multi-kaon (at least one kaon of the indicated charge 
and any other hadrons) production in neutrino and antineutrino scattering off iron and carbon. 
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tant. In other words, the antikaons seen in the detector 
are not the ones orginally produced in the first neutrino- 
nucleus interaction. With decreasing neutrino energy 
from high to moderate, the main source of secondary an- 
tikaons — the production of kaon-antikaon pairs in pion 
rescattering in nucleus — becomes more and more impor- 
tant in comparison with the initial production. For the 
MINOS neutrino flux, which peaks at about E v « 3 GeV, 



these secondaries just compensate the absorption of an- 
tikaons initially produced at higher energies. For the 
MINOS antineutrino flux the initial antikaon production 
is larger (because the flux has a larger part of higher- 
energy component), so the absorption dominates. For 
the NOVA flux, which peaks at about E v = 2 GeV, the 
cross section for antikaon production is an order of mag- 
nitude lower, but the effect of secondary production is 
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more pronounced. For kaons both the absorption and 
secondary production play a less significant role. 

The interplay between kaon absorption and production 
is illustrated in Fig. [TH which shows the kaon and an- 
tikaon cross sections at 3 GeV, i.e., the energy where the 
MINOS flux has its maximum. For K + the FSI show 
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FIG. 14. (Color online) Kinetic energy distributions per tar- 
get nucleon for K + and K~ production in neutrino scattering 
off iron at E v = 3 GeV. The various curves are explained in 
the figure; they show the interplay between absorption and 
secondary production of kaons and the influence of a hadron 
formation time. 

the expected behavior: the long-dashed line gives the 
spectrum of primordial kaons, i.e., those produced in the 
initial interaction, and the short-dashed lower curve gives 
the spectrum as it would look like if only absorptive in- 
teractions (in this case, only charge exchange into K° is 
possible because of strangeness conservation) took place; 
this is the result a typical Glauber or optical model cal- 
culation would give. The full result differs: Kaons with 
kinetic energies beyond about 0.4 GeV are suppressed by 
the FSI (solid curve) compared to the initial production 
distribution (long-dashed). This suppression is due to 
rescattering with accompanying energy loss; this explains 
the strong rise at lower kinetic energies. These results do 
not change when the formation time is set to because 



at this energy all produced hadrons (kaons and pions) 
are still slow enough to see the major part of the tar- 
get nucleus even when they start interacting only later. 
For K~ the situation is considerably more complicated: 
Absorption alone reduces the cross section significantly, 
and, due to the strong rise of the total K~N cross section 
at small momenta the cross section goes to zero at small 
momenta. The presence of secondary interactions in the 
FSI leads to a sizable enhancement of the yield at about 
0.2 GeV kinetic energy to a level well above the initial 
production. In this case the influence of formation time 
is still small, but visible; if the K~ start to interact ear- 
lier they will be more likely absorbed. Also contributing 
to this effect is that K~~ can be produced only from the 
sea and are thus nonleading hadrons so their interactions 
start with a lower cross section 12j. 



The Mineri^A experiment plans measurements of ex- 
clusive strangeness production. Both nuclear effects and 
experimental thresholds make it very difficult to distin- 
guish truly exclusive from semi-inclusive events and to 
isolate the true vN — > KX vertex. Here the same sec- 
ondary production can contribute to the cross section. 
To extract the genuine neutrino-strangeness cross sec- 
tions will require a very reliable modeling of the reaction 
mechanism. 

To summarize this section on particle production we 
give the integrated cross sections averaged over the MI- 
NOS and NOz/A fluxes in Table Q] 

VII. SUMMARY 

In this paper we have extended the GiBUU model, 
which was used to describe various neutrino reaction 
channels — quasielastic scattering, resonance produc- 
tion, and single-pion background — to higher energies 
where DIS becomes relevant. Since GiBUU provides a 
realistic treatment of nuclear effects in initial- and final- 
state interactions and incorporates now all the relevant 
reaction mechanisms it is a good tool to investigate the 
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TABLE I. Integrated cross sections per target nucleoli (in 
1(T 38 cm 2 ) for MINOS (Fe) and NO^A (C) experiments for 
various final states in neutrino induced reactions 

MINOS NOuA 



final state 


w/o FSI with FSI 


w/o FSI with FSI 


17T + 


0.76 


0.50 


0.65 


0.50 


l7T° 


0.30 


0.30 


0.22 


0.23 


l7T" 


0.0081 


0.14 


0.0024 


0.057 


1 P 


2.1 


0.43 


1.5 


0.62 


In 


0.82 


0.14 


0.45 


0.15 


1K+ 


0.15 


0.18 


0.056 


0.061 


1K° 


0.095 


0.13 


0.027 


0.035 


IK 


0.040 


0.038 


0.0068 


0.0069 


IK- 


0.039 


0.037 


0.0065 


0.0065 


multi 7r + 


2.2 


1.6 


1.3 


1.1 


multi 7T° 


1.5 


1.3 


0.75 


0.69 


multi 7T~ 


0.88 


1.0 


0.34 


0.42 


multi p 


2.1 


2.9 


1.5 


1.8 


multi n 


0.83 


2.5 


0.45 


1.2 



get, uncertainties in treating the initial-state interactions 
in the DIS channel lead to uncertainty of approximately 
2% in the inclusive cross section, which is compared with 
the recent experimental results on iron. A study of par- 
ticle knock-out events has shown that FSI, such as ab- 
sorption and rescattering of the outgoing particles inside 
the nucleus, are very important and lead to a significant 
modification of the particle spectra. Such spectra should 
be experimentally observable with the help of the fine- 
grained Minerva detector. Our results clearly indicate 
that, for the MINOS and NOvA fluxes, the spectra of 
the outgoing pions peak below a pion kinetic energy of 
100 MeV. For nucleons and kaons the FSI also signifi- 
cantly increase the cross section at low kinetic energies. 
For kaons, secondary production is seen to be very impor- 
tant and contributes a major part of the observed cross 
section. It will, therefore, be extremely difficult to ex- 
tract the strangeness producing vN cross sections from 
experiments with nuclear targets. 
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